NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


THESIS 


A  LINEAR  MANEUVERING  MODEL  FOR 
SIMULATION  OF  SLICE  HULLS 

by 

William  J.  Wolkerstorfer 
September,  1995 

Thesis  Advisor:  Fotis  A.  Papoulias 

Approved  for  public  release;  distribution  is  unlimited. 

19960221  034 


L 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No.  0704-0188 


i 


Public  reporting  burden  for  this  collection  of  infomiation  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instruction,  searching  existing 
data  sotirces,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or 
any  other  aspect  of  this  collection  of  infoimation,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Infonnation 

Operatioas  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction 

Project  (0704-0188)  Washington  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  Wan/c)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

September  1995  Master’s  Thesis 

4.  A  LINEAR  MANEUVERING  MODEL  FOR  SIMULATION  OF 

SLICE  HULLS 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S)  Wolkerstorfer,  William  J. 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey  CA  93943-5000 

8.  PERFORMING 

ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSORING/MONTTORING 
AGENCY  REPORT  NUMBER 

1 1 .  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the 
official  policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Govemmoit. 

12a.  DISTRIBUnON/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  luilimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (maximum  200  words) 

A  maneuvering  mathematical  model  is  generated  for  a  baseline  SLICE  ship  by  calculating 
the  linear  hydrodynamic  coefficients  using  existing  semi-empirical  methods  for  the  pods. 

Struts  are  modelled  as  flat  plates  attached  to  the  pods.  Pod/strut  combinations  are  translated  to 
the  ship's  reference  point  generating  die  ship's  linear  hydrodynamic  coefficients.  Eigenvalues 
of  the  horizontal  equations  of  motion  are  used  to  evaluate  the  stability  of  the  base  line 
configuration.  Rudder  effects  are  included  and  steering  characteristics  are  investigated  for 
various  rudder  aspect  ratios  and  locations.  Rudder  placement  is  discussed  in  the  context  of 
practical  considerations  and  model  test  results. 

14.  SUBJECT  TERMS  Ship  motion.  Equations  of  motion.  Simulators,  Hydrodynamics, 

Computerized  model,  hull  forms,  steering,  maneuvering 

15.  NUMBER  OF 

PAGES  74 

16.  PRICE  CODE 

17.  SECURTTY  CLASSIFI¬ 
CATION  OF  REPORT 
Unclassified 

18.  SECURITY  CLASSIFI¬ 
CATION  OF  THIS  PAGE 
Unclassified 

19.  SECURITY  CLASSFI- 
CATION  OF  ABSTRACT 
Unclassified 

20.  LIMITATION  OF 
ABSTRACT 

UL 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std.  239-1*  298-102 


1 


11 


Approved  for  public  release;  distribution  is  unlimited. 

A  LINEAR  MANEUVERING  MODEL  FOR  SIMULATION  OF  SLICE 

HULLS 


William  J.  Wolkerstorfer 
Lieutenant ,  United  States  Navy 
BME.,  University  of  Minnesota,  1988 

Submitted  in  partial  fiilfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 

from  the 


NAVAL  POSTGRADUATE  SCHOOL 
September,  1995 


Department  of  Mechanical  Engineering 


m 


IV 


ABSTRACT 


A  maneuvering  mathematical  model  is  generated  for  a  baseline  SLICE  ship  by 
calculating  the  linear  hydrodynamic  coefficients  using  existing  semi-empirical  methods  for 
the  pods.  Struts  are  modelled  as  flat  plates  attached  to  the  pods.  Pod/stmt  combinations 
are  translated  to  the  ship's  reference  point  generating  the  ship's  linear  hydrodynamic 
coefficients.  Eigenvalues  of  the  horizontal  equations  of  motion  are  used  to  evaluate  the 
stability  of  the  base  line  configuration.  Rudder  effects  are  included  and  steering 
characteristics  are  investigated  for  various  rudder  aspect  ratios  and  locations.  Rudder 
placement  is  discussed  in  the  context  of  practical  considerations  and  model  test  results. 


V 


TABLE  OF  CONTENTS 


L  INTRODUCTION . 1 

A.  MANEUVERING  AND  CONTROL . 2 

B.  THE  SLICE  HULL  FORM . 3 

H.  DETERMINATION  OF  HYDRODYNAMIC  COEFFICIENTS  OF  SLICE  PODS . 7 

A-  EQUATIONS  OF  MOTION  FOR  A  SUBMERGED  BODY  IN  THE  HORIZONTAL  PLANE: 
STEERING  SYSTEM  EQUATIONS _ 7 

B.  ESTIMATION  OF  HYDRODYNAMIC  COEFFICIENTS  FOR  A  BODY  OF  REVOLUTION 

USING  SEMI-EMPnaCAL  METHODS _ 11 

C.  SUBOFF  MODEL. _ 13 

D.  HYDRODYNAMIC  COEFnCIENT  CALCULATIONS  FOR  PODS _ 14 

E.  SLICE  POD  MODEL _ 19 

HL  HYDRODYNAMIC  COEFFICIENTS  OF  STRUTS _ 21 

A.  DESCRIPTION _ 21 

B.  CALCULATION  OF  THE  HYDRODYNAMIC  COEFFICIENT  CONTRIBUTIONS  FOR  A 

SINGLE  STRUT/POD  CONFIGURATION _ 22 

IV,  COMPLETE  HYDRODYNAMIC  DERIVATIVES  FOR  THE  SLICE  CONFIGURATION _ 25 

V.  RESULTS  AND  DISCUSSION _ 29 

A.  STEADY  STATE  TURNING  ABILITY _ 30 

B.  COURSE  KEEPING  AND  TURNING  ABILITY,  NOMOTO’S  FIRST  ORDER  MODEL _ 33 

VL  CONCLUSIONS . 37 

APPENDIX  A.  SAMPLE  CALCULATION  OF  THE  HYDRODYNAMIC  COEFFICIENTS  OF 
THE  SUBOFF  BODY  AND  SLICE  PODS. _ 49 

APPENDIX  B.  MATLAB  ROUTINE  FOR  CALCULATION  OF  THE  HYDRODYNAMIC 
COEFFICIENTS  FOR  A  SINGLE  STRUT/POD  CONFIGURATION . 55 

LIST  OF  REFERENCES _ 61 

INITIAL  DISTRIBUTION  LIST _ 63 


vii 


LIST  OF  FIGURES 


3 

1.  Conventional  SWATH  hull.  . . 

_  4 

2.  Profile  VIEW  of  SLICE  vessel.  [2] . 

3.  Stern  VIEW  OF  SLICE  [2] . ^ 

. 7 

4.  Body  fixed  axes  . . 

5.  Motion  stability  of  surface  ships.  [4] . 

6.  . . 

7.  Approximate  SUBOFF  body.  [10] . 

19 

8.  Profile  OF  SLICE  pod . 

9.  Baseline  strut/pod  configuration  for  SLICE  model . 21 

10.  Tactical  Diameter  (in  ship  lengths)  versus  rudder  percent  area  for  fixed  aspect  ratio 

rudders.  Rudders  are  attached  to  the  trailing  edge  of  the  after  strut.  Asterisks 

indicate  a  rudder  span  of  6  ft . 

1 1  Controls  fixed  eigenvalues  versus  rudder  percent  area  for  constant  aspect  ratios.  Rudder 

IS  FIXED  TO  TRAILING  EDGE  OF  AFT  STRUT . 

12.  D/L  versus  RUDDER  PERCENT  AREA  FOR  CONSTANT  ASPECT  RATIO  RUDDERS.  RUDDERS  FIXED  TO 

OF  AFT  POD . 

13.  Controls  fixed  eigenvalues  versus  rudder  percent  area  for  constant  aspect  ratio  rudders. 

Rudders  fdced  to  rear  of  aft  pod . ^  ^ 

14.  Limiting  D/L  versus  aspect  ratio  for  pod  and  strut  mounted  rudders . 42 

15.  Limiting  D/L  versus  rudder  percent  area  for  pod  and  strut  mounted  rudders . 43 

16.  Dominant  eigenvalue  versus  rudder  percent  area  for  constant  aspect  ratio  rudders  fixed 

TO  STRUT . 

17.  Dominant  eigenvalue  versus  rudder  percent  area  for  constant  aspect  ratio  rudders  fixed 

TO  AFT  POD . 

18.  D/L  VERSUS  PERCENT  AREA  FOR  RUDDER  WITH  DEADWOOD  ATTACHED  TO  POD . 46 

19.  1/K  VERSUS  1/T  FOR  CONSTANT  AREA  RUDDERS  ATTACHED  TO  TRAILING  EDGE  OF  AFT  STRUT . 47 

20.  l/K’-l/T’  DIAGRAM  FOR  CONSTANT  AREA  RUDDERS  FIXED  TO  TAIL  OF  AFT  POD . 48 


IX 


1.  INTRODUCTION 


The  SLICE  vessel  concept  was  bom  out  of  the  small  waterplane  area  twin  hull 
(SWATH)  hull  configuration  which  have  better  stability  characteristics  than  a  conventional 
mono  hull  of  similar  displacement  [1].  The  SWATH  hull  design  consists  of  two  submerged 
pods  which  provide  the  buoyancy  of  the  vessel  and  may  contain  various  pieces  of 
propulsion  or  auxiliary  machinery,  fuel  and  water  tanks  or  storage  spaces.  The  pods  may 
be  cylindrical  bodies  of  revolution,  similar  to  a  bare  submarine  hull,  or  may  be  some  other 
shape.  In  either  case  the  major  portion  of  the  vessel  static  stability  comes  from  the 
submerged  pods.  The  pods  are  connected  to  the  working  area  of  the  ship  by  means  of 
either  single  or  multiple  stmts.  At  this  point  in  time  there  are  no  standard  SWATH 
configurations.  Because  all  of  the  buoyancy  comes  from  the  pods  and  not  the  connecting 
stmts  they  must  have  a  large  lateral  separation  to  achieve  static  stability.  This  large 
separation  results  in  a  ship  with  a  much  wider  beam  and  shorter  length  than  that  of  a 
conventional  mono-hull  of  equal  displacement.  As  already  mentioned  the  stmts  provide 
very  little  additional  buoyancy  which  result  in  the  pods  being  deeply  submerged  giving  a 
much  deeper  draft  than  that  of  a  conventional  mono-hull  of  similar  displacement.  [1] 

The  primary  benefit  of  the  SWATH  hull  form  is  the  excellent  seakeeping  ability 
resulting  from  the  deep  draft  and  wide  beam.  The  SWATH  hull  form  has  pitch,  roll  and 
heave  characteristics  of  vessels  of  much  larger  displacements  providing  for  a  better  man- 
vessel  interface  in  heavier  seas.  The  man- vessel  interface  is  defined  as  the  ability  of  crew 
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members  to  perform  all  the  vessel  functions  without  suffering  from  the  effects  of  sea¬ 
sickness  or  being  swept  overboard  by  a  rogue  wave,  and  trying  to  maintain  one’s  balance 
while  working  in  heavy  seas.  The  wide  beam  may  leave  larger  open  deck  space  giving 
the  planner  or  designer  more  flexibility  in  placement  of  deck  fittings  or  machinery.  The 
vessel’s  shorter  length  may  give  the  captain  a  better  overall  view  of  the  working  deck  area 
allowing  swifter  and  more  efficient  command  of  operations  on  deck.  While  the  shorter 
length  may  be  considered  as  an  advantage  it  comes  at  a  cost  of  tendency  to  pitching  due 
to  a  relatively  small  longitudinal  metacentric  height  (GMl).[1] 

A.  MANEUVERING  AND  CONTROL 

In  addition  to  its  superior  seakeeping  ability,  the  SWATH  hull  is  also  directionally 
very  stable,  a  direct  result  of  the  deep  drafts  and  large  underwater  hull  profile  area. 
Conventional  mono-hulls  may  require  the  addition  of  a  rudder,  fin  or  deadwood  in  order 
to  achieve  directional  stability  whereas  the  SWATH  concept  hulls,  which  are  in  essence 
large  fins  or  wings,  are  exceptionally  stable  directionally  at  low  speeds  and  increasingly 
so  as  flow  over  the  hull  increases.  This  increased  stability  also  has  its  drawback  because 
along  with  the  ability  to  keep  a  straight  course  in  seas  or  wind,  it  requires  a  much  larger 
yawing  moment  to  turn.  At  lower  speeds  this  large  yaw  moment  can  be  overcome  with 
maneuvering  thrusters  or  by  differential  thrast  between  the  engines.  At  higher  speeds  the 
differential  thrust  is  not  as  effective  and  the  rudders  must  provide  the  necessary  moment 
to  initiate  the  turn.  [1] 
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Rudders  in  use  in  current  SWATH  concept  ships  are  standard  balanced,  horn  or 
spade  rudders  and  are  placed  aft  as  in  a  conventional  mono-hull  design.  They  are 
attached  to  the  strut  or  haunch  above  and  behind  the  screw  as  shown  in  Figure  1.  The 
rudder’s  effectiveness  is  increased  by  placing  it  in  the  propeller  wash.  [1] 


Figure  1.  Conventional  SWATH  hull.  [1] 

B.  THE  SLICE  HULL  FORM 

The  SLICE  hull  differs  from  the  SWATH  hull  in  that  the  SLICE  consists  of  two 
pods  on  each  side  providing  buoyancy  and  stability  versus  one  as  in  the  SWATH  hulls. 
Figure  2  shows  a  profile  view  of  the  SLICE  with  major  dimensions  and  without  rudders. 
Figure  3  is  a  stem  view  ,  diagram  of  the  SLICE  vessel  without  radders  or  stabilizing  fins, 
showing  fore  and  aft  pod  offsets.  In  terms  of  maneuvering  ability  the  SLICE  presents 
some  challenges  in  the  location  of  the  rudders.  Should  the  mdders  be  placed  in  the 
propeller  wash  increasing  the  effectiveness  of  the  radder  but  reducing  the  turning 
moments  created  by  such  a  rudder?  Should  they  be  placed  on  the  after  pod?  The  after 


3 


strut?  Perhaps  on  the  leading  edge  of  the  aft  strut  to  place  it  somewhat  in  the  propeller 
wash  but  increasing  the  moment  it  can  create  or  even  on  the  leading  edge  of  the  fore  strut. 


41* 


Figure  3.  Stem  view  of  SLICE  [2]. 
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II.  DETERMINATION  OF  HYDRODYNAMIC  COEFFICIENTS  OF 

SLICE  PODS 

A.  EQUATIONS  OF  MOTION  FOR  A  SUBMERGED  BODY  IN  THE 

HORIZONTAL  PLANE.  STEERING  SYSTEM  EQUATIONS 

The  equations  of  motion  for  a  submerged  body  in  the  horizontal  plane  are  similar 
in  nature  to  that  of  an  airplane  in  flight.  The  vehicle  is  assumed  to  be  in  steady  state 
motion  in  calm  seas  without  the  effects  of  currents  and  waves.  All  steering  system 
dynamics  terms  are  referenced  to  a  body  fixed,  orthogonal  coordinate  system  with  the 
origin  fixed  at  the  ship’s  reference  point,  which  is  usually  fixed  at  the  amidships  main 
deck  but  can  be  chosen  to  be  at  any  point.  The  x  axis  is  fixed  to  the  longitudinal  center  of 
the  ship  and  is  positive  forward.  The  y  or  lateral  axis  is  positive  to  starboard  and  the  z 
axis  is  positive  in  the  downward  direction  relative  to  the  ship.  [3,4,5] 


Figure  4.  Body  fixed  axes  [5] 

For  a  vehicle  with  a  high  aspect  ratio  such  as  a  surface  ship,  the  forces  are 
predominantly  forces  of  lift  developed  as  a  result  of  the  relative  motion  of  the  vessel  and 
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the  fluid.  The  longitudinal  forces  of  drag,  while  important  for  the  determination  of 
power  requirements  are  neglected  for  steering  system  dynamics.  It  is  assumed  that  no 
loss  of  speed  takes  place  in  maneuvering.  In  straight  line  motion  the  lift  forces  on  each 
side  of  a  body  which  is  symmetric  about  a  longitudinal  axis  are  equal  and  sum  to  zero. 
When  the  fluid  impinging  on  the  vessel  has  an  angle  of  attack  with  respect  to  the 
longitudinal  axis,  then  the  lift  forces  are  no  longer  in  balance  and  dynamic  forces  are 
developed.  In  other  words  when  the  vessel  has  a  side  slip  velocity  component  (v)  in  the 
y-direction  an  angle  of  attack  relative  to  the  fin  or  strut  is  caused.  The  angle  of  attack  (a) 
is  defined  as  v/U  for  small  v  and  it  gives  rise  to  the  lift  forces  and  moments  on  the  vessel 
hull.  The  ability  of  the  hull  to  maintain  a  straight  line  path,  in  other  words  straight  line 
stability,  will  be  measure  a  of  the  stability  of  that  particular  hull  form.  The  linearized 
sway  and  yaw  equations  of  motion  for  a  vessel  in  the  horizontal  plane  are  fully  developed 
in  Principles  of  Naval  Architecture  [3,4]  and  repeated  here  in  matrix  form  as:  [3] 


m-F;  -Yr  ' 

m 

V 

1 

1 

■s. 

j 

V 

_-Ni  L-N;. 

• 

.r. 

LiVv  Nr  J 

r 

Equation  1  is  often  more  conveniently  represented  in  matrix  form  as: 

Mx  =  Ax  (2) 
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The  components  of  the  M  and  A  matrices  of  equation  2  represent  the  linearized 
hydrodynamic  derivatives  of  the  forces  developed  by  sideslip  (v)  and  yawing  (r)  motions 
of  the  ship.  Coefficients  m’  and  Izz’  are  the  non-dimensionalized  mass  and  mass  moment 
of  inertia  of  the  vessel  about  the  z  axis.  The  forces  developed,  in  their  linearized  form  are 
represented  as  the  derivative  times  the  motion.  In  other  words  Yy’v  represents  the  lateral 
force  developed  by  the  side  slip  velocity  v  and  Yr’r  represents  the  lateral  force  developed 
by  the  yawing  motion  r.  In  similar  fashion  the  moments  Ny’v  and  Nr’r  represent  the 
moments  on  the  vessel  about  the  selected  reference  point  resulting  from  the  lateral 
velocity  and  yawing  motion.  [3] 

In  order  to  fully  study  the  steering  system  dynamics  the  effects  of  the  rudder  must 
be  added  to  equation  2  as  follows. 

Mx  =  Ax+B6  (3) 

Where  B  is  a  vector  containing  the  linearized  lateral  force  and  moment 
contributions  of  the  rudders  as  shown  in  equation  3  and  6  is  the  angular  displacement  of 
the  rudder  (-  for  right  radder,  +  for  left  rudder  for  standard  stem  radder). 

B  =  N'a]  (4) 

In  determining  the  controls  fixed  stability  of  a  hull  shape  it  is  the  eigenvalues  of 
equation  2,  re-written  as 
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that  represent  a  measure  of  the  controls  fixed  stability  of  the  ship.  Controls  fixed  stability 
is  the  stability  associated  with  a  straight  line  path  of  the  vessel  with  zero  rudder  angles. 
Complex  conjugate  eigenvalues  with  negative  real  parts  correspond  to  a  directionally 
stable  vessel  but  with  oscillatory  response  characteristics  as  shown  in  Figure  5  case  HA. 
Real,  negative  eigenvalues  will  be  directionally  stable  with  no  oscillations  as  in  case  IIB. 
One  or  more  positive  eigenvalues  indicate  a  hull  form  that  is  directionally  unstable.  [3,4] 


CASE  I 

5rNAIGHT>I.INE  STAWim’ 


CASE  HA 

EMAECTIONAL  stahlitv 
LOW  DAAWING 


CASE  IIB 

OIAECTIONAL  STABILITY 
high  DAMTING 


CASE  III 

KISITIONAL  STABILITY 


Figure  5.  Motion  stability  of  surface  ships.  [4] 


B.  ESTIMATION  OF  HYDRODYNAMIC  COEFFICIENTS  FOR  A  BODY  OF 
REVOLUTION  USING  SEMI-EMPIRICAL  METHODS 

There  are  several  different  methods  for  estimation  of  the  hydrodynamic 
coefficients  for  a  submerged  bodies  of  revolution  relying  on  some  semi-empirical  relation 
'  to  account  for  viscous  or  vortex  effects  on  the  body.  The  methods  considered  for  the 

SLICE  pods  are  summarized  in  a  concise  and  comprehensive  report  from  the  Naval 
Coastal  System  Center  [6].  It  includes  methods  from  United  States  Air  Force  Data 
Compendium,  (USAF  DATCOM),  Naval  Coastal  Systems  Center  (NCSC)  and  Naval 
Ship  Research  and  Development  Center  (NSRDC).  While  each  agency  has  published 
several  reports  on  the  subject,  Peterson  [6]  outlines  the  strengths  and  weaknesses  of  each 
of  the  above  mentioned  methods.  Comparison  of  the  various  methods  is  made  for  both 
submarine  shaped  and  torpedo  shaped  bodies  of  revolution  developed  by  each  source  [6]. 
The  acceleration  hydrodynamic  coefficients  will  be  calculated  by  the  methods  of 
Humphreys  and  Watkinson  [7]. 

All  of  the  methods  discussed  are  easily  calculated  both  by  hand  and  are  amenable 
to  computerized  calculations  [6].  However  the  methods  of  DATCOM  and  NCSC  rely 
primarily  on  geometric  considerations.  Methods  that  required  body  masses  or  mass 
distributions  such  as  NSRDC  methods  are  discounted  at  this  preliminary  point  in  the 
modeling  process  as  they  are  not  as  well  defined  or  determined  as  the  physical 
dimensions  of  the  SLICE  pods.  [6] 
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Because  of  a  lack  of  captive  model  tests  or  full  scale  maneuvering  trial  data  it  is 
necessary  to  build  the  computer  model  of  the  SLICE  hull  form  in  several  parts.  The  first 
part  is  the  pod  itself  which  resembles  a  submarine  form  and  is  fully  described  below. 
Attached  to  each  of  the  pods  is  the  strut  connecting  the  pod  to  the  box.  Struts  are 
modeled  as  a  flat  plate  fixed  fin  attached  to  the  pod.  The  hydrodynamic  coefficients  of 
each  strut/pod  combination  will  be  calculated  individually  for  each  pod  with  the  reference 
point  for  moment  coefficients  taken  at  the  pod  mid  length  point.  The  effects  of  each  of 
the  four  strut/pod  configurations  will  then  be  non-dimensionalized  with  respect  to  the 
vessel  length  between  perpendiculars  (Lbp)  and  translated  to  the  ships  reference  point  and 
combined  to  get  the  total  coefficients  for  the  ship.  [3] 

While  the  SLICE  pod  form  is  not  exactly  that  of  a  submarine  the  semi-empirical 
methods  used  to  calculate  the  coefficients  of  the  pods  will  first  be  verified  with  captive 
model  test  data  for  the  SUBOFF  body  [8,9]  by  idealizing  the  body  model.  The  idealized 
model  body  will  consist  of  three  simple  geometric  shapes  and  is  discussed  below  with  the 
SUBOFF  body. 

The  NCSC  method  is  based  on  data  presented  in  Reference  [10]  and  gives  results 
consistent  with  the  DATCOM  calculations  [6]  for  the  SUBOFF  body  [8, 9].  The 
DATCOM  method  is  used  for  slice  pod  calculations  because  of  ease  of  programming  for 
parametric  studies  without  requiring  coefficients  and  correction  factors  derived  from 
graphical  techniques.  Results  of  the  calculations  presented  in  [6]  for  methods  not  used 
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but  considered  will  be  presented  for  the  SUBOFF  [8]  body  without  outlining  the 
calculation  procedures  or  details. 


C.  SUBOFF  MODEL 

David  Taylor  Research  Center  (DTRC)  has  designed  the  Defense  Advanced 
Research  Projects  Agency  (DARPA)  SUBOFF  model  5470  and  5471  described  in 
Reference  [9].  The  hydrodynamic  coefficients,  as  determined  by  captive-model  tests  are 
shown  in  [8]  and  will  be  used  for  eomparison  purposes.  In  this  comparison  the  referenee 
point  (Xm)  is  located  at  the  body’s  center  of  gravity  as  opposed  to  it’s  amidships  length  as 
used  for  the  slice  pods.  The  planform  of  the  DARPA  model  5470  is  shown  in  Figure  6 
[10]. 


Figure  6.  SUBOFF  Body  [9] 

NCSC  [10]  uses  a  simplified  body  of  revolution  for  estimation  of  the 
hydrodynamic  coefficients  which  very  closely  resembles  the  SLICE  pod  planform.  It 
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consists  of  an  elliptical  nose,  a  cylindrical  or  parallel  mid  body  section  and  a  conical  base 
section  [10].  Figure  7  shows  the  simplified  body  used  for  estimation  of  the 
hydrodynamic  coefficients. 


Figure  7.  Approximate  SUBOFF  body.  [10] 


D.  HYDRODYNAMIC  COEFFICIENT  CALCULATIONS  FOR  PODS 

Y  y  :  Calculation  of  the  body  alone  normal  force  coefficient  using  the  DATCOM 

and  NCSC  methods  are  similar  in  nature  and  form  but  vary  in  the  method  whereby  the 
axial  position  on  the  body  where  the  flow  becomes  predominantly  viscous  is  measured. 
The  Air  Force  DATCOM  method  is  based  on  aircraft  and  missile  data  and  is  measured 
form  the  nose  of  the  body  and  is  determined  from  the  empirical  relation.  [6] 

l.  =  0m,  +  0£21l_  (6) 

where  L  is  the  distance  from  the  nose  of  the  vessel  to  the  point  where  the  flow  becomes 
predominantly  viscous,  1b  is  the  overall  length  of  the  body  and  l^s  is  the  distance  from  the 
nose  of  the  body  to  the  point  of  maximum  slope  on  the  tail  section.  In  the  case  of  a 
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conical  tail  the  point  of  maximum  slope  is  at  the  tip  making  l^s  =  1b  in  this  case  and  for 
the  SLICE  pod.  [6] 


The  lift-curve  slope  is  computed  from 

=  (7) 

Ub 

where  Ci„  is  the  lift  slope  curve,  k2  and  ki  are  the  Lamb’s  coefficients  of  inertia  for  a 
prolate  ellipsoid  in  axial  and  cross  flow  respectively.  Sy  and  Sb  are  the  cross  sectional 
areas  of  the  base  at  point  L  and  Sb  is  the  maximum  cross  sectional  area  of  the  body.  [6] 
The  normal  force/angle  of  attack  curve  slope  is  calculated  as 

(yvL=^[c,.+cj  w 


where  Cdo  is  the  drag  coefficient  of  an  elliptical  form  at  a  zero  angle  of  attack  and  1  is  the 
total  length  of  the  pod.  [6] 

^ :  Body  alone  yawing  moment  coefficient  calculations  presented  by  Peterson 

[6]  show  very  good  correlation  with  experimental  data  for  all  methods.  The  NCSC  and 
NSRDC  calculations  gave  the  best  correlation  with  the  SUBOFF  body.  However,  the 
data  presented  in  [9]  is  too  limited  for  SLICE  pods  and  the  NSRDC  method  [6]  requires 
knowledge  of  the  mass  of  the  body.  The  DATCOM  method  is  used  as  follows  [6]. 


(9) 
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where  Cma  is  the  yawing  moment  curve  slope  and  is  calculated  as 


where  S  is  the  cross  sectional  area  of  the  body  as  a  function  of  x  (S=S(x))  and  Xm  is  the 
moment  reference  point  for  the  submerged  body  [6]. 


y  ,  ^  :  Body  Alone  rotary  derivatives  calculated  with  the  DATCOM  and 


NCSC  methods  are  identical  in  form  except  for  the  calculation  of  the  point  of  viscous 
flow  separation  as  discussed  above.  The  DATCOM  method  is  [6]. 


where  V  is  the  displacement  volume  of  the  submerged  body  and  Ic  is  the  distance  from 
the  nose  to  the  center  of  buoyancy.  Stb  is  the  area  of  the  truncated  base  ly .  [6] 

Acceleration  hydrodynamic  coefficients  used  for  the  SLICE  pods  are  based  on  the  works 
of  Humphreys  and  Watkinson  [7]  and  calculated  as  follows 
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(13) 


(14) 


where  kb  is  Lamb’s  coefficient  of  added  moment  of  inertia  for  a  prolate  ellipsoid  body. 
Izdf  is  the  mass  moment  of  inertia  of  the  displace  fluid  about  the  z  axis.  Terms  Icb  and  leg 
are  distance  from  the  nose  to  the  center  of  buoyancy  and  center  of  gravity  [7].  Table  1 
shows  the  results  of  the  above  calculations  for  the  SUBOFF  body. 
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Experiment 


DATCOM/%  error 


NCSC  /  %  error 


NSRDC/%  error 


-0.005948 

-.005836 

-2.475% 

-0.007106 

19.47% 

-.007555 

27.02% 

N\ 

-0.012795 

-.013550 

5.90% 

0.012693 

-199.2% 

.012693 

-199.2% 

Yr 

0.001811 

-0.001473 

■181.34% 

-0.002038 

-212.53% 

-0.003259 

-279.96% 

Nr 

-0.001597 

-0.001183 

-25.92% 

-0.001628 

1.94% 

-0.001405 

-12.02% 

Y\  * 

-0.013278 

-0.015272 

15.02% 

n/a 

n/a 

Y\ 

.000060 

0 

0 

0 

N- 

.000202 

0 

0 

0 

N'r 

-0.000676 

-0.000694 

2.66% 

n/a 

n/a 

Table  1.  Comparison  of  semi-empirical  methods  of  determining  hydrodynamic 

coefficients. 
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E.  SLICE  POD  MODEL 

The  model  used  to  represent  the  SLICE  pods  is  closely  approximated  by  a  36  ft 
long  8  ft  diameter  body  of  rotation  in  the  format  described  above.  It  consists  of  an 
elliptical  nose  6  ft  long,  with  a  parallel  mid-body  section  14  ft  long  with  a  conical  base  16 
ft  long  as  shown  in  Figure  8. 


Calculations  of  the  hydrodynamic  coefficients  of  the  SLICE  pods  are  done  using 
methods  of  [6]  and  [7]  for  each  pod  without  fins  and  based  on  a  reference  point  (Xm)  at 
the  mid  length  of  the  pod.  Reference  point  x^  for  the  pods  is  eighteen  feet  (1b/2).  Pod 
and  strut  coefficients,  when  calculated  individually  will  be  non-dimensionalized  with 
respect  to  pod  length  (1b)  and  later  re  non-dimensionalized  with  respect  to  the  vessel 
length  between  perpendiculars  (Lbp). 
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III.  HYDRODYNAMIC  COEFFICIENTS  OF  STRUTS 


A.  DESCRIPTION 

The  struts  connecting  the  pods  to  the  body  of  the  SLICE  vessel  are  modeled  as  flat 
plates  twenty  four  feet  in  length,  six  feet  in  depth  with  a  zero  taper  ratio  and  zero  sweep 
back  angle  (A).  Figure  9  shows  the  baseline  strut/pod  configuration  used  in  the  modeling 
process. 


Figure  9.  Baseline  stmt/pod  configuration  for  SLICE  model. 

The  hydrodynamic  derivatives  for  the  struts  are  determined  in  two  parts  for  each 
strut.  The  first  part  is  the  fin  section  formed  where  the  base  begins  to  taper  and  the 
second  is  the  six  by  twenty-four  portion  which  makes  up  the  majority  of  the  strut.  Each 
part  is  modeled  as  a  flat  plate  fixed  fin  located  at  a  distance  Xf  from  the  center  of  the  pod 
to  the  quarter  chord  point  of  the  fixed  fin.  Once  the  hydrodynamic  coefficients  are 
calculated  for  each  strat/pod  configuration  individually  they  are  translated  to  the  reference 
point  of  the  vessel  and  added.  At  this  point  in  the  modeling  process  the  coefficients  for 
the  pod  and  strut  are  non-dimensionalized  with  respect  to  the  length  of  the  pod.  [3] 


B.  CALCULATION  OF  THE  HYDRODYNAMIC  COEFFICIENT 

CONTRIBUTIONS  FOR  A  SINGLE  STRUT/POD  CONFIGURATION 

The  contributions  to  the  hydrodynamic  derivatives  of  a  fixed  fin  are  derived  in 
[3,4]  and  result  in  the  following 


Y.f  =  -A 


\ 


+  C.J 


(15) 


Where  A’  is  the  non-dimensional  area  of  the  fin,  dCi/da  is  the  lift-angle  of  attack  (a) 
slope  curve.  Cd  is  the  profile  drag  coefficient  for  a  flat  plate  at  zero  angle  of  attack  and  is 
taken  as  zero  for  fins  and  rudders.  Subscript  f  denotes  a  fin  or  rudders.  For  large  ships 
with  relatively  small  rudders  or  fins  the  location  of  the  point  of  application  of  the  forces  is 
ignored  and  Xf  is  taken  as  the  distance  from  the  ship’s  reference  center  to  the  half  chord 
point  of  the  fin  [3,4].  Because  of  the  large  size  of  the  struts  relative  to  the  length  of  the 
vessel  it  is  felt  that  this  additional  distance  cannot  be  ignored.  Theoretically,  the  point  of 
application  of  the  hydrodynamic  force  on  a  flat  plate  fin  with  zero  taper  angle  is  at  the 
quarter  chord  point  [1  l].so  that  when  taking  the  moment  contribution  of  the  fixed  fin 
from  the  center  of  the  pod  Xf  is  the  distance  from  Xm  to  the  quarter  chord  point .  The 
distance  x’f  is  the  non-dimensional  distance  from  the  pod  reference  point  to  the  quarter 
chord  point  i.e.  Xf/L.  The  contributions  of  the  fixed  fin  to  the  pod  coefficients  are  given 
in  equations  3  through  5  and  are  easily  derived  by  summing  forces  and  moments  about 
point  Xm.  [3] 


N’vf  =  Y’vfX’f 

(16) 

Y\f  =  x'fY\f 

(17) 

N\  =  Y\x^f 

(18) 

The  acceleration  coefficient  of  a  flat  plate  under  lateral  acceleration  is 
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Coefficients  a  and  ao  in  equations  2  and  6  are  the  effective  aspect  ratio  and 
effective  geometric  aspect  ratio  respectively.  For  a  fixed  fin  which  is  in  contact  with  a 
hull  or  groundboard  such  that  flow  is  not  allowed  around  the  tip  or  root  chord  the 
effective  aspect  ratio  is  twice  the  aspect  ratio  for  the  fm.  In  the  case  of  the  SLICE  struts 
which  pass  through  the  water  surface  and  are  in  contact  with  the  pod  the  effective  aspect 
ratio  is  taken  as  2a.  [3,4] 

Following  similar  logic  for  the  added  mass  or  acceleration  derivatives  are  shown 
in  equations  7  through  9. 


mf=T;fx\  (20) 

yv=rv/jcv  (21) 

=  (22) 

The  hydrodynamic  coefficients  contributions  of  the  strut  can  now  be  added 
algebraically  to  that  of  the  pod  as  follows.  Subscripts  p  and  f  designate  the  coefficient 
contributions  from  the  pod  and  fin  or  strut  respectively.  [3] 


Y’v  =  Y’v)p  +  Y’vf 

(23) 

Y’r  =  Y’,)p  +  Y’rf 

(24) 

N’v  =  N’v)p  +  N’vf 

(25) 

N’,  =  N’,)p  +  N’rf 

(26) 

Similar  logic  is  followed  for  the  added  mass  coefficients  and  is  repeated  for  each  of  the 
four  strut/pod  configurations.  [3] 
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IV.  COMPLETE  HYDRODYNAMIC  DERIVATIVES  FOR  THE 

SLICE  CONFIGURATION 

Translating  the  hydrodynamic  coefficient  contributions  of  each  of  the  four 
strut/pod  configurations  is  performed  by  taking  force  and  moments  about  the  reference 
point  of  the  vessel.  Principles  of  Naval  Architecture  [3]  derives  the  equations  with 
reference  to  an  arbitrary  point  on  a  conventional  hull  and  they  are  adapted  here  for  a 
single  strut/pod  as  follows. 

Y\  =  Y\  (27) 

rr  =  Y\x\  +  Y\  (28) 

where  x’f  here  represents  the  distance  from  the  amidships  reference  point  to  the  reference 
point  of  the  pods  with  the  correct  signs  (+  for  forward  pods,  -  for  after  pods) 

N'.  =  r.x',-^N',  (29) 

N>r,x',+ry,*N\x',^N,  (30) 

The  acceleration  derivatives  are  determined  with  similiar  logic  and  are  not  shown 

here. 

The  final  step  in  determining  the  overall  hydrodynamic  derivatives  is  to  non- 
dimensionalize  with  respect  to  the  characteristic  length  of  the  vessel  (Ibp)  and 
algebraically  sum  all  of  the  contributing  components  together  [3].  Assuming  that  the 


25 


forces  on  each  of  the  fore  and  aft  strut/pods  are  equal  port  and  starboard  we  have  for  the 
baseline  configuration  the  coefficients  as  follows. 


(Fv)  ={(Fv)^+(yv)J 

(Afv)=2[(Ar,),,„+(7Vv)J 

(iV.)=2[(Ar.)^+(7Vr)J 


Similiar  procedure  is  followed  for  the  added  mass  derivatives  and  will  not  be  shown  here. 

The  results  of  the  above  calculations  are  summarized  in  Table  2.  Note  that  the 
non  dimensionalization  of  the  pod  alone  and  the  two  strut/pod  configurations  are  with 
respect  to  the  length  of  the  pod.  The  base  line  vessel  is  non  dimensionalized  with  respect 
to  Lbp- 
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SLICE  Pod 

Alone 

Fwd  strut/pod 

configuration 

Aft  Strut/pod 

configuration 

Base  line  vessel 

w/o  rudders. 

Y\ 

-0.016244 

-0.019149 

-0.020316 

-0.078930 

Tr 

-0.002498 

-0.002518 

-0.001135 

-0.004044 

N\ 

-0.038754 

-0.005614 

-0.004231 

-0.016428 

N\ 

-0.006396 

-0.000576 

-0.000257 

-0.010332 

Ti 

-0.044189 

-0.012682 

-0.012982 

-0.051328 

T; 

0.0 

-0.001346 

0.005860 

0.005617 

m 

0.0 

-0.001305 

-0.005860 

-0.001945 

-0.001799 

-0.000221 

-0.000070 

-0.00564 

Table  2.  Summary  of  hydrodynamic  coefficients  for  SLICE  vessel. 
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V.  RESULTS  AND  DISCUSSION 


After  the  hydrodynamic  coefficients  for  the  base  line  configuration  are  calculated 
the  stability  eigenvalues  of  the  equations  of  motion  are  calculated  to  be  negative  real 
values  -0.089  and  -1.208.  These  values  indicate  that  the  base  line  configuration  without 
fins  is  inherently  stable  due  to  the  large  fin  area  associated  with  the  struts  in  agreement 
with  the  SWATH  results  [1].  A  parametric  study  was  then  conducted  to  study  the  effects 
of  moving  the  struts  toward  and  away  from  the  reference  point.  The  results  were  as 
expected,  the  eigenvalues  become  more  negative  indicating  increased  vessel  stability  and 
a  related  increase  in  the  moment  necessary  to  initiate  a  turn.  [1,3] 

Three  rudder  configurations  are  studied.  The  first  is  a  pair  of  rudders  attached  to 
the  trailing  edge  of  the  after  strut  without  deadwood  or  attaching  structures.  The  second 
is  a  rudder  attached  to  the  rear  of  the  after  pod  without  any  deadwood  and  the  third  is  a 
rudder  as  a  portion  of  a  deadwood  attached  to  the  rear  of  the  after  pod.  At  this  point  no 
attempt  is  made  to  simulate  any  interference  effects  from  fins,  fixtures,  hinges  or  any 
additional  structure  necessary  to  affix  the  rudder  or  deadwood  to  the  vessel.  Additionally 
no  effort  is  made  to  establish  the  pressure  effects  between  the  pods  and  their  effect  on  the 
rudders.  The  size  and  aspect  ratio  of  the  rudders  is  varied  for  each  configuration  studied 
with  the  limiting  rudder  geometry  being  the  span.  The  maximum  span  is  limited  to  6  feet 
because  it  is  felt  that  this  would  be  the  largest  span  rudder  that  could  be  put  in  place  on 
the  aft  strut.  The  aft  strut  has  a  clearance  of  8.5  feet  between  the  water  surface  and  the 
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pod.  This  model  also  assumes  that  the  ship  does  not  pitch  or  roll  thereby  keeping  the 
entire  mdder  submerged. 

A.  STEADY  STATE  TURNING  ABILITY 

Once  the  inherent  stability  of  SLICE  is  determined,  the  effects  of  mdders  are 
introduced.  The  rudders,  modeled  as  flat  plates,  are  not  affixed  to  any  structure  at  the 
root  and  tip  allowing  flow  across  the  tip  and  root  of  the  mdder.  Each  mdder  is  assumed 
to  have  a  full  effect  on  the  turning  ability  because  of  the  large  separation  between  aft 
stmts  and  neglecting  pressure  effects  between  them.  Calculation  of  the  hydrodynamic 
coefficients  of  the  mdder  and  the  stmts  are  performed  by  the  same  procedure  with  two 
notable  exceptions.  First,  the  reference  point  for  mdder  effect  is  measured  from  and 
calculated  with  respect  to  the  vessel  reference  point.  This  effectively  translates  the  force 
and  moment  effects  of  the  mdder  directly  to  the  reference  point  of  the  ship.  The  second 
exception  is  based  on  the  assumption  that  the  mdder  does  not  broach  the  surface  or  come 
in  close  contact  with  any  stmcture  at  the  tip  or  root  chord.  This  makes  the  effective 
aspect  ratio  simply  a  vice  2a  as  discussed  above  [3]. 

Figure  lOshows  the  tactical  diameter  in  ship  lengths  (D/L)  versus  percent  mdder 
area  for  several  constant  aspect  ratio  mdders  attached  to  the  trailing  edge  of  the  aft  stmt. 
The  percent  underwater  (u/w)  profile  area  is  the  percentage  of  mdder  area  (Ar)  with 
respect  to  design  length  between  perpendiculars  and  draft  (LbpT)  or  A/(LbpT).  The  aspect 
ratios  in  Figure  10  increase  as  curves  go  down.  The  top  curve  is  for  an  aspect  ratio  of  1 
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and  the  lower  curve  is  an  aspect  ratio  of  5.  The  asterisks  (*)  represent  the  limiting  rudder 
span  of  6  feet  and  increases  for  points  moving  right  on  the  curves.  Rudders  with  a  low 
aspect  ratio  will  tend  to  be  short  and  fat  extending  further  aft  while  a  higher  aspect  ratio 
rudder  will  tend  to  be  tall  and  slender  as  shown  in  Figure  10. 

The  results  indicate  that  this  model  will  have  tactical  diameters  between  4  and  7 
ship  lengths  depending  on  rudder  area  and  aspect  ratio.  Conventional  mono  hulls  and 
current  SWATH  configurations  have  rudder  areas  of  between  1  and  3  percent  of  u/w 
profile  area  [1].  For  rudder  areas  in  this  range  Figure  10  indicates  that  aspect  ratios 
between  1  and  3  would  be  viable  alternatives  without  consideration  for  rudder  stock 
moments,  rudder  bending  stresses  or  rudder  protrusions  aft  of  the  vessel  stem. 

Figure  1 1  shows  the  controls  fixed  eigenvalues  versus  percent  rudder  area  for 
several  constant  aspect  ratio  rudders  attached  to  the  aft  strut.  In  each  case  the  aspect  ratio 
is  increasing  as  the  curves  go  down.  The  upper  most  curve  is  for  an  aspect  ratio  of  1  and 
the  lower  most  curve  for  an  aspect  ratio  of  5.  Because  the  eigenvalues  increase  (become 
more  negative)  with  increasing  mdder  size  the  vessel’s  straight  line  stability  is  increasing 
at  the  cost  of  decreased  turning  ability. 

In  the  case  of  a  rudder  attached  to  the  rear  of  the  aft  pod  the  resulting  steady  state 
turning  radius  and  eigenvalue  trends  are  similar  to  those  discussed  above.  The  only 
difference  is  in  the  magnitude  of  the  values.  As  expected  a  radder  on  the  pod,  located 
further  from  the  ships  center  of  gravity  will  have  an  increased  turning  moment  and 
thereby  decrease  the  turning  diameter  of  the  ship.  Figure  12  shows  that  the  mdder  on  the 
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aft  pod  will  generate  a  smaller  turning  circle  due  to  the  larger  rudder  moment  and  Figure 
13  shows  the  controls  fixed  eigenvalues  for  the  second  rudder  configuration.  Figure  14 
and  Figure  15  compares  the  minimum  tactical  diameters  for  the  limiting  cases  represented 
by  the  asterisks  of  figures  10  and  12.  The  results  shown  in  Figure  10  through  Figure  13 
are  reasonably  good  qualitative  descriptions  of  the  stability  but  do  not  clearly  indicate  the 
interaction  between  the  turning  ability  and  course  stability  of  SLICE. 

Biancardi  [12]  in  a  study  of  a  data  base  of  173  vessels  of  various  configurations 
quantified  the  dominant  eigenvalue  for  vessels  which  displayed  satisfactory  turning 
ability[12].  It  shows  that  the  dominant  eigenvalue  for  both  the  bare  hull  and  hull  and 
rudder  configurations  were  stable  and  had  good  turning  ability  if  the  dominant  eigenvalue 
was  between  -.331  and  zero.  Dominant  eigenvalue  less  than  -.331  indicates  stable  hulls 
but  with  reduced  turning  ability.  Vessels  with  positive  eigenvalues  were  unstable 

[3,4,11].  [12] 

Figure  16  and  Figure  17  show  the  dominant  eigenvalue  (Oa)  versus  rudder  percent 
area  for  strut  mounted  and  pod  mounted  radders  respectively.  The  asterisk’s  show  the 
limiting  span  as  described  above.  Based  on  [12]  the  curves  indicate  that  both  rudder 
configurations  will  result  in  adequate  stability  and  turning  ability  for  all  aspect  ratios  up 
to  and  beyond  the  limiting  span  selected. 

Case  three,  a  rudder  with  deadwood  attached  to  the  pod  is  predicted  based  on  the 
turning  characteristics  in  the  previously  discussed  figures.  The  deadwood  is  modeled  as  a 
flat  plate  fixed  fin  at  the  aft  pod  with  an  eight  foot  span  corresponding  to  the  maximum 
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diameter  of  the  pod  and  an  aspect  ratio  of  three  yielding  a  deadwood  chord  of  2.67  feet. 
Rudder  size  for  this  configuration  is  limited  to  a  span  equal  to  that  of  the  deadwood. 
Figure  18  is  the  D/L  versus  percent  rudder  area  for  the  third  case  and  clearly  show  that  for 
a  deadwood  attached  to  the  aft  pod  the  turning  ability  is  markedly  decreased.  The 
dominant  eigenvalue  for  this  configuration  is  -.321  for  all  aspect  ratio’s  and  rudder  areas 
indicating  that  the  vessel  is  stable  but  will  have  poorer  turning  performance  when 
compared  the  first  two  configurations  with  greater  (less  negative)  eigenvalues. 

B.  COURSE  KEEPING  AND  TURNING  ABILITY,  NOMOTO’S  FIRST 
ORDER  MODEL 


Nomoto  presented  a  simplified  analysis  of  steering  system  dynamics  by 
manipulating  the  linearized  steering  system  equations  of  motion.  He  developed  the  first 
order  equation  for  the  yaw  rate,  given  in  terms  consistent  with  modem  control 
engineering  given  in  equation  35.  [4] 


where 


T= 


K'= 


r.-N; 

N\ 

Ks 

N\ 


(36) 

(37) 
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The  indices  K’  and  T’  represent  the  ratios  of  non-dimensional  coefficients  as  the  yaw 
inertia  /yaw  damping  coefficient  and  the  turning  moment/yaw  damping  coefficients 
respectively.  [4] 

Equation  35  shows  that  the  yaw  rate  r’  can  be  represented  as  a  first  order  system 
with  gain  K’  and  time  constant  T’ .  In  steady  state  the  yaw  rate  per  length  (r’  =  r/L)  is 
simply  K’5.  A  larger  gain,  K’  provides  a  greater  steady  state  turning  ability  i.e.  increased 
steady  state  yaw  rate,  while  an  increased  value  of  1/T’  (smaller  time  constant)  will 
provide  a  quicker  response  to  the  rudder.  In  addition  to  indicating  the  response  to  the 
rudder  T’  can  also  provide  an  indication  of  the  ships  course  stability.  Since  a  quick 
response  to  rudder  order  is  desired  in  a  course  keeping  /  control  situation  an  increased 
value  of  T’  is  also  a  quantitative  indication  of  the  ship’s  course  stability.  [4] 

Figure  19  shows  the  turning  system  gain  (1/K’)  and  time  constant  (1/T’)  for 
constant  rudder  areas  based  on  the  Nomoto  first  order  approximation  for  a  rudder  on  the 
aft  strut.  Each  curve  represents  a  constant  area  rudder  with  the  higher,  steeper  curves  to 
the  left  representing  small  area  while  the  low,  shallow  curves  to  the  right  are  higher 
rudder  area.  The  area  limits  represented  are  between  five  square  feet  and  thirty  five 
square  feet  in  five  square  foot  increments  with  the  maximum  span  of  six  feet.  Rudder 
span  decreases  as  the  points  move  upward  on  the  curve  with  a  corresponding  increase  in 
rudder  chord  and  decreasing  aspect  ratio.  The  data  in  Figure  19  show  that  for  a  given 
area  rudder,  a  decrease  in  aspect  ratio  as  we  move  up  the  curve  will  decrease  the  turning 
rate  (decrease  in  K’)  of  SLICE  but  with  a  corresponding  increase  in  the  sensitivity  to 


34 


rudder  commands  and  course  stability.  It  is  also  readily  seen  that  increasing  the  rudder 
area  will  increase  both  course  stability  and  responsiveness  to  rudder  commands. 
However  in  terms  of  rudder  response  and  course  stability  the  range  of  1/T’  shown  is  very 
small  indicating  that  the  course  stability  is  relatively  constant  for  a  large  range  of  rudder 
sizes  and  corresponds  well  with  the  eigenvalue  view  of  stability  discussed  previously. 

Figure  20  is  the  first  order  gain  and  time  constant  (1/K’  -1/T’)  curve  for  case  2. 
The  primary  difference  between  the  curves  in  Figure  19  and  Figure  20  is  found  in  the 
value  of  T  for  a  given  rudder  area.  A  rudder  attached  to  the  pod  will  have  smaller  time 
constant  than  a  rudder  attached  to  the  strut  indicating  a  better  response  to  rudder 
commands.  Additionally  this  shows  that  for  either  rudder  configuration  of  a  given  area, 
course  stability  is  not  very  sensitive  to  changes  in  rudder  size. 
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VI.  CONCLUSIONS 


Based  on  the  linearized  model  presented  herein  the  SLICE  hull  has  a  very  good 
turning  ability  as  well  as  course  stability  as  shown  by  the  steering  system  eigenvalues  and 
first  order  Nomoto  model.  The  results  discussed  clearly  indicate  that  the  turning  ability  is 
enhanced  with  a  rudder  place  further  aft  or  made  larger  but  as  in  all  engineering 
endeavors,  trade  offs  and  practical  considerations  must  be  taken  in  to  consideration.  Low 
aspect  ratio  rudders  (aspect  ratio  1  or  2)  should  be  considered  only  with  caution  as  they 
may  extend  too  far  aft,  potentially  beyond  the  stem  of  the  ship.  They  will  require  more 
power  to  operate  because  of  the  larger  moments  about  the  mdder  stock  and  will  have 
much  larger  bending  stresses  associated  with  them.  Very  high  aspect  ratio  mdders 
(aspect  ratios  four  to  five  or  greater)  should  also  be  considered  carefully  as  they  may  often 
broach  the  surface  of  the  water  in  heavy  seas,  decreasing  their  effectiveness.  The  data, 
based  on  the  model  presented  supports  the  use  of  a  stmt  mounted  mdder  of  aspect  ratio 
between  two  and  four. 

As  a  recommendation  for  furthur  research,  we  should  analyze  full  scale  sea  trail 
data  of  the  SLICE  hull  and  identify  the  hydrodynamic  coefficients.  This  will  enable  us  to 
compare  the  linearized  analytical  model  with  experimental  results. 
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Diameter/Lenglh 


Percent  Area 


;ure  10.  Tactical  Diameter  (in  ship  lengths)  versus  rudder  percent  area  for  fixe 
io  rudders.  Rudders  are  attaehed  to  the  trailing  edge  of  the  after  strut.  Asteris 
licate  a  rudder  scan  of  6  ft. 


APPENDIX  A.  SAMPLE  CALCULATION  OF  THE 
HYDRODYNAMIC  COEFFICIENTS  OF  THE  SUBOFF  BODY  AND 
SLICE  PODS. 
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Rotary  Coefficients  using  DATCOM  method 

Planform  sectional  area  of 

S  ijj  :  =  2-  r(x)  dx  truncated  base  (a.k.a.  effective  base  in  NCSC 

Jl  terminology) 
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Acceleration  Hydrodynamic  Coefficients  of  pod/body 
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Normal  force  Coefficient  using  Datcom  method 


C  Do  ~  ^-29  from  Hoemer  page  1  - 1 8  for 
hemispherical  heads  and  noses 


^ms  ^ms  14.292*ft 
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point  where  viscous 
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area  of  effective  base 
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Nyprime  Coefficient  using  DATCOM  method 
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determination  of  center  of  buoyancy 

Assuming 

-  Pod  has  weight  uniforaily  distributed  over  length 


x-7i-2  r(x)  ydx 


7i-r(x)  •2*ydx 


l^b=6.492-ft 


this  is  assuming  that  the  body  is 
homogenous  solid 


1  cb  “  6.59003 '  ft  length  to  center  of  buoyancy  as  given  by 

SUBOFF  document  [8,9] 


Munk  coefficients 


A 

16^4  n 


\  2-e^  /  \  1  -  e 


3  2  1-e 


2-Po 


k2  =0.950  kj  =0.026 


2-  •  2-e^-  2  - 


(Po-«o) 


kb  =0.854 


Other  area  and  mass  calculations 


m  :  =  .018078  from  SUBOFF  data 


S(x)  : = 71- r(  cross  sectional  area  of  body 
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d  :=1.667  fl 
lpMB-7.3125.ft 

Ij,  :=  3.6459  ft 

Ip 

Xm  :  =  6.604- ft 


V  :  =  1 .059- 1 0"  ^-ft^-sec'  ^  y  ;  =  62.4  — 

Length  of  ellipsoidal  nose  assumed  to  end  at  front  of  sail 
Length  of  parallel  mid  body 

Length  of  base  assumed  a  conical  shape 

length  between  perpendiculars  1  g  =  1 4.29 1 7  *11 

length  of  forebody  1  p  =  10.646*ft 

position  of  reference  center  for  calculations 


7 

g 


^  bp  ^  B 


numpts  -  1 000  number  of  points  used  for  plotting  and  such 


Body  geometric  (radius)  characteristics  and  cross  sectional  areas 
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V=23.462*ft^ 


max  cross  sectional  area  of  body 


0-ft 
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Sp  =  19.593 -ft^ 


cross  sectional 
Planform  area 
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APPENDIX  B.  MATLAB  ROUTINE  FOR  CALCULATION  OF  THE 
HYDRODYNAMIC  COEFFICIENTS  FOR  A  SINGLE  STRUT/POD 
CONFIGURATION. 


function  [HCP]=hcsp(L,c,b,xf,fa) 

% 

%  Hydrodynaniic_Coefficients_Slice_Pod 

% 

%  .m  File  to  calculate  the  hydrodynamic  coefficients  of  the  SLICE 
%  pods  based  on  geometric  considerations  using  semi-empirical 
%  methods  developed  by  various  sources,  (refer  to  main  main 
%  document  for  descriptions  of  methods.  All  coefficients  are 
%  calculated  in  dimensionless  form  using  SNAME  convention. 

% 

%  Reference  documents  are  Naval  Coastal  Systems  Center  report  (NCSC 
%  TM-291-80)  "Evaluation  of  Semi-Empirical  Methods  for  Predicting 

%  Linear  Static  and  Rotary  Hydrodynamic  Coefficients"  Authored  by 

%  R.  S.  Peterson. 

% 

%  Reference  for  acceleration  Coefficients  is 

%  Naval  Coastal  Systems  Lab  "Prediction  of  Acceleration 
%  Hydrodynamic  Coefficients  for  Underwater  Vehicles  from  Geometric 
%  Parameters"  (NCSL  TR-327-78)  by  D.E.  Hymphreys  and  K.W.  Watkinson 
% 

%  Input  Arguments 
% 

%  L  Column  vector  of  lengths  defining  body  of  rotation 

%  c,b  Chord  and  span  of  fiat  plate  strut. 

%  xf,xa  Length(ft)  from  reference  point  to  1/4  chord  of  strut 
%  fa  integer  flag  (1  or  2)  to  indicate  wether  fore  or  aft 

%  strut/pod  combination  is  being  calculated 

% 

%  Output  Arguments 


% 

% 

HCP  Vector  of  values  as  follows 

% 

HCP(l) 

Yv 

% 

HCP(2) 

Yr 

% 

HCP(3) 

Nv 

% 

HCP(4) 

Nr 

% 

HCP(5) 

Yvdot 

% 

HCP(6) 

Yrdot 
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%  HCP(7)  Nvdot 

%  HCP(8)  Nrdot 

%  Coefficients  Yrdot  and  Nvdot  are  assumed  zero  for  pods 
%  because  they  are  generally  very  small  and  contribute  little. 

% 

%  **********  Start  of  Executable  code  ********** 

% 

load  slice 

% 

%  Fill  output  vector  for  pod  alone  non-dimensionalized  with  respect  to 
%  length  of  the  pod  (Calculated  in  mathcad  document  "fwdpod.mcd"  and 
%  "aftpod.mcd"  which  include  the  effects  of  the  tail  as  a  flat  plate. 

% 

if(fa==  1) 

% 

%  Foreward  pod  coefficients  with  tail  fin  as  flat  plate  calculated 
%  independantly. 

% 

Yv=-.017390;  Yr=-.002371;  Nv=-.038626;  Nr=-.006268; 

Yvdot=-.044430;  Yrdot=0.000027; 

Nvdot=.000127;  Nrdot=-.001802; 

elseif  (fa  ==  2) 

% 

%  Aft  pod  coefficients  with  tail  fin  as  flat  plate 

% 

Yv=-.023406;  Yr=-.001105;  Nv=-.037361;  Nr=-.005003; 

Yvdot=-.047953;  Yrdot=0.000732; 

Nvdot=.001393;  Nrdot=-.001941; 

end 

% 

HCP(l)=Yv; 

HCP(2)=Yr; 

HCP(3)=Nv; 

HCP(4)=Nr; 

HCP(5)=Yvdot; 

HCP(6)=0; 

HCP(7)=0; 

HCP(8)=Nrdot; 

% 

%  Calculate  the  coefficients  for  strut  as  fixed  fin  attatched  to  pod 

% 

HCS=hcfps(xf,c,b,T,LB); 

% 
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%  combine  and  non-dimensionalize  the 
%  between  perpendiculars. 

% 

HCP(1)=HCP(1)*(LB/Lbp)'^2 

HCP(2)=HCP(2)*(LB/Lbp)'^3 

HCP(3)=HCP(3)*(LB/Lbp)^3 

HCP(4)=HCP(4)*(LB/Lbp)M 

HCP(5)=HCP(5)*(LB/Lbp)^3 

HCP(6)=HCP(6)*(LB/Lbp)M 

HCP(7)=HCP(7)*(LB/Lbp)M 

HCP(8)=HCP(8)*(LB/Lbp)'^5 


coefficients  wit  vessel  length 


+  HCS(1)*(LB*T)/Lbp'^2; 

+  HCS(2)*(LB''2*T)/Lbp''3; 

+  HCS(2)*(LB''2*T)/Lbp''3; 

+  HCS(4)*(LB'^3*T)/LbpM; 

+  HCS(5)*(LB'^2*T)/Lbp''3; 

+  HCS(6)*(LB'^3*T)/LbpM; 

+  HCS(7)*(LB''3*T)/LbpM; 

+  HCS(8)*(LBM*T)/Lbp'^5; 
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function  [HCS]=hcfjps(xf,c,b,T,L) 


% 

%  Hydrodynamic_Coefficients_Flat_Plate_Strut 

% 

%  .m  file  to  calculate  the  hydrodynamic  coefficients  of  the  SLICE 
%  vessel  struts  as  flat  plates  placed  at  a  position  of  xf  forward  or 
%  aft  of  the  pod's  reference  point.  This  routine  views  the  struts  as 
%  a  fixed  fin  attatched  to  the  pod  with  the  pressure  center  acting  at 
%  the  point  Cmc4  fore  (+)  or  aft  (-)  of  the  1/4  chord  point  of  the  fin. 

% 

%  Coefficients  are  non-dimensionalized  as  for  a  surface  vessel  iaw 
%  SNAME  conventions. 

% 

%  INPUT  ARGUMENTS 

% 

%  xf  Location  of  center  of  strut  relative  to  vessel's  center 
%  reference  point. 

%  c  Average  chord  of  strut/flat  plate  (width) 

%  b  Span  of  flat  plate/strut 

%  T  Draft  of  vessel  (for  non-dimensionalization) 

%  L  Length  of  vessel  (for  Non-dimensionalization) 

%  OUTPOUT  ARGUMENTS 
% 

%  HCS  Vector  of  8  hydrodynamic  coefficients  in  the  following  order 

%  HCS(l)  Yv 

%  HCS(2)  Yr 

%  HCS(3)  Nv 

%  HCS(4)  Nr 

%  HCS(5)  Yvdot 

%  HCS(6)  Yrdot 

%  HCS(7)  Nvdot 

%  HCS(8)  Nrdot 

% 

%  LOCAL  VARIABLES 


% 

% 

a 

Aspect  ratio 

% 

ag 

Geometric  aspect  ratio 

% 

A 

Area 

% 

Cmc4 

Quarter  chord  moment  coefficient 

% 

Xf 

Point  of  reference  for  moments 

% 

%  **********  Start  of  Executable  Code  ********** 
% 
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%  Initialize  variables 
% 

A=b*c;  %  Area  of  strut 

ag=2*(b/c);  %  Aspect  Ratio  of  fin 

a=ag;  %  geometric  aspect  ratio 

Cla=(0.9*pi*2*a)/(sqrt(a'^2+4)+1.8);  %  Simplified  version  of  eq  23a  (pna) 
Cmcl=l/2-((  1.11  *sqrt(ag''2+4)+2)/(4*ag+2)) 

Cmc4=(.25-Cmcl)*(Cla) 

Xf=xf; 

xf,Xf,Cmc4 

Yv=-abs(A*Cla/(L*T)); 

Yr=(Xf/L)*Yv; 

Nv=Yr; 

Nr=(Xf/L)''2*Yv; 

Yvdot=-(2*pi*c*b^2)/((L'^2*T)*sqrt(ag'^2+l)); 

Yrdot=-(Xf/L)*Yvdot; 

Nvdot=(Xf/L)*Yvdot; 

Nrdot=(Xf/L)''2*Yvdot; 

% 

%  Fill  output  vector 
% 

HCS(l)=Yv; 

HCS(2)=Yr; 

HCS(3)=Nv; 

HCS(4)=Nr; 

HCS(5)=Yvdot; 

HCS(6)=Yrdot; 

HCS(7)=Nvdot; 

HCS(8)=Nrdot; 
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